The sizes of galaxy halos in galaxy cluster Abell 1689 
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ABSTRACT 

The multiple images observed in galaxy cluster Abell 1689 provide strong constraints not only 
on the mass distribution of the cluster but also on the ensemble properties of the cluster galaxies. 
Using parametric strong lensing models for the cluster, and by assuming well motivated scaling 
laws between the truncation radius s and the velocity dispersion a of a cluster galaxy we are able 
to derive sizes of the dark matter halos of cluster galaxies. 

For the scaling law expected for galaxies in the cluster environment (s oc a), we obtain 
s = 64l}^ x [a J 220 km/s) kpc. For the scaling law used for galaxies in the field with s oc a 2 



-14 

we find s = 66^ 



16 



[a I 220 km/s) 2 kpc. Compared to halos of field galaxies, the cluster galaxy 



halos in Abell 1689 are strongly truncated. 

Subject headings: gravitational lensing - cosmology:dark matter - galaxies:clusters:individual:Abell 1689 



1. Introduction 

Although galaxies are the units of objects seen 
on cosmological distances, very little is known 
observationally about the extent of dark matter 
halos surrounding the galaxies beyond the light 
emitted by the gas and stars in them. Based on 
numerical simulations these dark matter halos are 
expect ed to extend out to several hundred kpc 
(e.g. Tormen et al. 1998). Rotation curves of 



19891 : IPrada et al.ll2003h . This has recently been 
done for the Milky Way bv lBattaglia et alJ (|2005h 



who were able to measure the radial velocity dis- 
persion profile of the Galaxy out to 120 kpc. The 
method works only for field galaxies since it is sen- 
sitive to other nearby massive galaxies and hence 
the g alaxies studied need to be isolated (jBrainerd 
200411 . 



spiral galaxies and the line of sight velocity dis- 
persions of the stars in elliptical galaxies can be 
mea sured only out to some tens of kiloparsecs (see 
e.g. ISofue fc Rubii]|200lt iBender et al]|l994l and 
references therein). The results from these two 
methods indicate that the masses of galaxies con- 
tinue to grow roughly linearly with the radius, i.e. 
the matter in galaxies is closely isothermal. The 
radial velocities of satellites of galaxies can also be 
used to estimate the masses of the i r host galaxies 
(e.g. Hartwick k. Sargent] 1978 : Zaritsky et al 



Gravitational lensing is an ideal tool to measure 
the extents of dark matter halos around galaxies 
since no optical tracers within the halo are needed. 
Instead, the halo mass can be inferred from the 
gravitational lensing of background sources. 

Weak lensing can be used to study galactic 
dark matter halos statistical ly. The field starte d 
from the pioneering work of Tvson et al.l (1984), 
and galaxy-galaxy lensing h as now been success- 



fully used both in the field (Braiucr d et al 



dell' Antonio fc Tvsonl 1 19961: iHudson et al 



1996; 



1998; 



1 



Fischer et al.l2000t | Smith et al.ll200lHWilson et al 



2001 



2004(1 and in clusters (INataraian et alJ 199 



McKay et all200lllHoekstral2003l:lHqekstra et al 



Geiger & Schneider 1999; Nataraian et al. 
Gavazzi et al. 20041 Limousin et al.l 20061 ) to 




2002 
mea- 
sure the masses and extents of galaxy halos. The 
signal is very weak for individual galaxies and 
needs to be collected from many galaxies, pos- 
sibly adopting various scaling laws to compare 
measurements from lensing galaxies with different 
luminosities. The different works generally find a 
tangential shear 7 that decreases like 7 cx 1/9 with 
the radius 9, i.e. the halos stay roughly isother- 
mal beyond the luminous component. In the field 
the signal from galaxies has been measured out to 
~200kpc(e.g. IWilson et alj200ltlHoekstra et al 
20041) . 



The galaxy truncation in clusters has been 
stu died both theoretic a lly a nd observationally 



INataraian fc Kneibl (Il997l); INataraian et al 



.19981; Geiger fc Schneider! (ll998lll999h:lNataraiari et al 
( 20021) ; bavazzi et al.l (|20ol ; lLimousin et al.l(|2005T ~ 
20061). Strong t r uncati o n of galaxies is found in 



Nataraian et alJ (Il998l). INataraian et all (|2002h 
and Limousin et al.l ( 20061 ) when compared to 
galaxies in the field (truncation radii s* of L* 
galaxies span 17-55 kpc in the 6 clusters vs 
s*=2 64±42 kpc//i7 in the field, iHoekstra et al 



2004). There is a general agreement that the ha- 
los of galaxies in dense environments are truncated 
relative to those in the field although the uncer- 
tainties are still large and the sample of clusters 
used is inhomogeneous. The inherently statisti- 
cal nature of the methods and the need to as- 
sume certain scaling laws further complicate the 



case. The method used in INataraian et al. (|l998i 



2002) also requires that the parameters of the 
smooth cluster componen t are known accurately 
(jNatarajan fc Kneiblll997l ) . This is achieved by in- 
corporating also s trong lensing con straints in the 
clusters enabling iNatarajan et al. to accurately 
model the cluster profile. 

The typical radius of an Einstein ring of a 
galaxy is at most a few arcseconds and so in the 
strong lensing regime it is not possible to probe 
the extent of dark matter halos beyond a few arc- 
seconds directly. This makes strong lensing unfea- 
sible to study the dark halos of individual galaxies 
beyond a few arcseconds in the field. In clusters 
of galaxies however the combined potentials of the 



smooth dark matter halo of the cluster as a whole 
and those of the individual galaxy halos are re- 
sponsible for the lensing of background sources. 
This enables us to statistically probe galaxy halos 
in dense environments using strong gravitational 
lensing. Since lensing constrains only the total 
potential of the cluster it is important to have a 
large number of multiple images with a large radial 
spread over the cluster in order to investigate the 
different mass components separately. Abell 1689 
is ideally suited for this task with the large num- 
ber of identified multiple images of many different 
b ackgro und sources and well defined strong lensing 



models (Broadhurst et al. 2005; Diego et al. 200 



Zekser et al.l 120061 lHalkola et all 120061 ). In this 



p aper we use the stron g lensing models developed 



Halkola et alJ ( 20061 ) to study the truncation of 
cluster galaxy dark matter halos in A1689. We 
demonstrate that the models are indeed sensitive 
to the total mass contained in the cluster galaxies 
and derive sizes for the galaxies in the cluster. 
This is the first time the sizes of galaxy halos have 
been measured in dense cluster environments with 
strong lensing only. 

In section [21 we give a br i ef sum mary of the 
models used in lHalkola et al.l (|2006) in particular 
the modeling of the galaxy component of the clus- 
ter, in section [3] we outline the methodology used 
to study the truncation of the cluster galaxies. 
The results are presented in section 2] and in sec- 
tion [5] we perform several checks to demonstrate 
that the results obtained are robust and reason- 
able. In section [6] we compare the results with 
earlier published studies of galaxy halo truncation 
before concluding in section [7J 

The cosmology used throughout this paper is 
ft m =0.30, O\=0.70 and H o =70 km/s/Mpc. 

2. Strong Gravitational Lensing Model for 
A1689 

The strong lensing models in thi s work are 
based on the parametric models used in lHalkola et al.l 
(2006) to study the mass profile of A1689 in detail. 
Here we give a short summary of t he strong lens- 



i ng m odeling but refer the reader to lHalkola et al 
(|2006l) for the details 



The multiple images in Halkola et al. ( 20061 ) 



2 



were i n most part those identified in lBroadhurst et al 
(|2005j ). In total a 107 multiple images in 31 multi- 
ple image systems and 1 arc were used. In 5 cases 
at least one of the images in a system had also a 
spectroscopic redshift and the redshifts were kept 
fixed for these systems. The redshifts of another 
26 multiple image systems were estimated using 
photometric redshifts. In these cases the redshift 
of an image system was allowed to find its best 
redshift within the estimated photometric redshift 
errors. The arc was too faint for good photometry 
and its reshift is left unconstrained. 

The mass in the cluster is assumed to be in two 
smooth DM halos that are described by either non- 
singular isothermal ellipsoids (NSIE) or elliptical 
Navarro- Frenk- White profiles (ENFW) . The small 
scale mass structures associated with the galax- 



ies ar e modeled with BBS profiles (jBrainerd et al 
1996f ). The BBS profile is a singular isothermal 



sphere with a truncation radius s. In the central 
regions (r < s) the density profile is isothermal 
(p oc r~ 2 ) but there is a truncation of the halo at 
radius s after which the density falls sharply with 

-4 

p oc r . 

The velocity dispersions of the galaxies are es- 
timated using the Fundamental Plane (FP). The 
FP tics kinematic (velocity dispersion), photo- 
metric (effective surface brightness) and morpho- 
logical (half light radius) galaxy properties to- 
gethe r ( Dressier et all 1987 ; Diorgovski fc Davisl 



19871 : iBender et al.l I1992T ). Measuring morpho 



logical and photometric properties of the galax- 
ies allows us to estimate the galaxy kinematics. 
We assume that the central velocity dispersions of 
galaxies, as derived from the FP, are equal to the 
halo velocity dispersions, and that the masses in 
disks can be neglected. For some fainter galax- 
i es we have also used th e Faber- Jackson relation 
( Faber k. Jackson 19761 here after FJ relation) 
that relates the absolute magnitude of a galaxy 
to its velocity dispersion. The FJ relation has a 
large intrinsic scatter and the velocity dispersion 
obtained using FJ have a larger uncertainty than 
the ones obtained with the FP. 

The truncation radii of the galaxies are as- 
sumed to follow a scaling law of the form s ga i = 
{vgai/ '(J°) a • In this paper we dis c uss th e 



s° x 



same scaling laws as in lHalkola et al.l ([2006), 
namely a = 1 and a = 2. a— I corresponds 
to tidal truncation of halos in dense cluster en- 



vironment (jMerrittl 119831 ) whereas galaxies with 
a=2 have a constant mass-to-light ratio and 
is us ually assumed in weak lensing analyses 



Brainerd et al 



19961: iNataraian et al.|[l99l 



Hoekstra et al.l 120041 ). In this paper we explore 
further the radial exte nt of the galaxy halo s for 
the scaling laws used in lHalkola et al. ( 20061 ). 

A1689 is an excellent candidate for this work 
since the large number of multiple images ensures 
not only that the global mass profile can be con- 
strained very accurately but also the relative con- 
tributions of the smooth DM and galaxies can be 
determined as will be shown later. 

3. Methodology 

In this paper take advantage of the unique op- 
portunity presented in A1689 to use strong lens- 
ing and the significant contribution of the cluster 
galaxies on the positions of the impressive number 
of multiple images observed in the cluster. The ef- 
fect is only observable in the total fit quality and 
is hence statistical in nature in that extension of 
individual galaxies cannot be determined. 

Unlike the usual galaxy-galaxy lensing in which 
foreground galaxies weakly distort the shapes of 
background galaxies this method relies on the 
changes induced by galaxies on the positions of 
multiply imaged background galaxies. This strong 
galaxy-galaxy lensing is only applicable in the 
strong lensing regime where multiple images are 
observed over a large range of cluster centric radii 
so that they pose strong constraints both on the 
total cluster potential and also on the galaxies. 

In this section we outline the method used to 
measure the extents of galaxy halos. The strong 
lensing models are constrained by the observed 
multiple images. The positions of the images can 
be measured to an accuracy of better than 1 pixel 
or 0.05" on the images from the Advanced Cam- 
era for Surveys. The only other measurables are 
the redshifts of the cluster and the multiple im- 
ages. The redshift of the cluster is well estab- 
lished from spectroscop i c surveys (iTeague et al 



11990 : iBalogh et alJl2002t iDuc et alJl2002f ) and the 
overall mass scale of the cluster is fixed by the 5 
spectroscopic redshifts of multiple image systems. 
The major uncertainty in the models is the inclu- 
sion of the cluster galaxies. In the following we 
describe the Monte-Carlo simulations used to find 



3 



the normalization of the scaling law, s°, and how 
these simulations can also be used to estimate the 
error in s° due to the uncertainties in the observ- 
ables. 

3.1. Monte-Carlo Simulations 

A Monte-Carlo run consists of reassigning a new 
velocity dispersion, <tmc, to each cluster galaxy 
based on the value <T ga i and estimated error de- 
termined using the FP or the FJ relation. The 
new <tmc of a galaxy was drawn from a Gaussian 
distribution centered on a ga i with a width corre- 
sponding to the estimated error. The multiple im- 
age positions were similarly varied with assumed 
error of 1 pixel. In this way we have constructed a 
simulated galaxy that has properties similar to the 
one observed within our estimates of the errors. 

This cluster can now be analyzed in the same 
way as the 'original' cluster. This means that we 
find the optimal parameters for the two smooth 
DM halos (positions, ellipticitics, position angles 
and the two free parameters of the halos: veloc- 
ity dispersion and core radius for the NSIE profile 
and concentration and virial radius for the ENFW 
profile) and redshifts for the image systems with 
photometric redshifts. 

The simulated clusters u sed in this work are 
the same that were used in lHalkola et all (|200fih 
to derive errors for the total mass profile and the 
parameters of smooth DM halo. 

In this work we concentrate on the normaliza- 
tion of the truncation radius s° for two scaling 



laws, a 



Halkola et al 



1 a nd a = 
(|2006l) . 



2, which was not done in 
This means that in addi- 
tion to optimizing the above mention parameters 
we also find the optimal value of s° for each sim- 
ulated cluster. This is explained below. 

3.2. Determining s° for a Monte-Carlo 
Run 

The optimal s° for each Monte-Carlo run was 
taken as the one with the minimum ( x 2 ) 1 ^ 2 when 
s° was progressively increased from 20 kpc to 
200 kpc. The parameter optimization was per- 
formed by a source plane minimization for com- 
putational reasons. In all subsequent analysis we 
have used an image plane x 2 defined as the sum 
of the squared distances between the observed im- 
ages and ones predicted by our models. ( x 2 ) 1 ^ 2 



is hence the rms distance between the observed 
images and the corresponding model image posi- 
tions. 

We use the shapes of these ( x 2 ) 1 ^ 2 vs - s ° curves 
to convince the reader that there is clear signal and 
that s° can be constrained in clusters using strong 
lensing once sufficiently many multiple images can 
be used to constrain the models. 

3.3. Determining s° for the Cluster 

The shape and spread of the ( x 2 ) 1 ^ 2 vs - s ° 
curves could in principle also be used to derive con- 
fidence limits on s°. This, however, would require 
us to perform more simulations to derive appro- 
priate Ax 2 levels for the confidence limits. The 
( X 2 } 1 ^ 2 vs - s ° curves do however demonstrate 
that there is a strong and clean signal that can be 
used to derive s° and the errors for a given scaling 
law. 

The best fitting s° and the errors for the cluster 
are derived from the distribution of the s° values 
obtained in the Monte-Carlo runs instead. 

4. Results 

For the scaling law we need to choose a refer- 
ence cr°. The derived s° is then the truncation 
radius of a galaxy with a velocity dispersion equal 
to this a. The truncation radii of galaxies with 
different <rs can then be obtained using the appro- 
priate scaling law. In this work we simply assume 
a fiducial value of ct°=220 km/s. To compare the 
s° obtained in this work with literature one should 
scale our s° by ( of it / 220 km/s ) a . 

In creating the simulated clusters the veloc- 
ity dispersions for the galaxies are drawn from a 
Gaussian distribution and hence we do not expect 
to see significant differences in the shapes of the 
individual ( x 2 ) 1 ^ 2 vs - s ° curves between the dif- 
ferent Monte-Carlo runs that would arise from a 
systematic change in the galaxy component. The 
curves do vary in their absolute ( x 2 } 1 ^ 2 level how- 
ever. For this reason we have normalized the indi- 
vidual curves to their respective median ( x 2 ) 1 ^ 2 
in order to bring all the curves to a similar ( x 2 ) 1 ^ 2 
level. After this the curves have been renormal- 
ized to the level of the mean median ( x 2 ) 1 ^ 2 of all 
the curves. The scaling of the individual curves is 
necessary in order to combine the information on 
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s° from the different curves. 

In Fig. Q] the mean curves for 1000 simulated 
clusters for each of the smooth DM profiles used 
are shown. NSIE is shown as a dotted line and 
squares, ENFW as a dashed line and triangles. 
The left panel shows the curves for a = 1 (s oc a) 
and the right for a = 2 (s oc a 2 ). Combining the 
two smooth DM profiles yields the solid line and 
points shown as circles. The points show the 2- 
sigma clipped mean for each s°. The error bars 
show the dispersions of the final clipped points for 
a given s . 

The smaller scatter in the points for NSIE mod- 
els is an indication that the renormalized curves 
are very similar while the considerable scatter for 
the ENFW models shows that the curves differ 
not only in the absolute x 2 level but also in their 
shape. The combined curve has been calculated 
from the NSIE and ENFW curves and not from 
the curves of each individual Monte-Carlo run for 
the two models. ENFW smooth DM profiles gen- 
erally favor slightly larger values for s° than the 
NSIE models (As - 10 kpc). 

The { x 2 ) 1 ^ 2 vs - s ° curves are flatter at large 
(>50 kpc) s° than at smaller s°, this is especially 
apparent for the cases where the smooth DM is 
described by an NSIE profile. A possible reason 
for the shallower slope on the logarithmic hori- 
zontal scale (linear in fractional change in mass) 
in Fig. [1] is that the larger extent of the halos, 
and hence a smoother combined mass profile of 
the galaxies, makes it easier for the smooth DM 
component to compensate for the change in the 
mass in the galaxies. In the small s° regime the 
galaxies have significant local contribution to the 
image positions which cannot be easily compen- 
sated by the smooth DM component. 

The curves for the two scaling laws are very sim- 
ilar and we are not able to differentiate between 
them in terms of quality of fit. This is also seen in 
weak lensing determinations of the extensions of 
the dark matter halos of galaxies (Limousin 2006, 
private communication). A possible explanation 
is that instead of measuring the extension directly 
we are in fact measuring the mass of the galaxies. 
This then creates a degeneracy between the two 
parameters in the scaling law, namely a and s°. 



For a different value of a the same total mass in 
the galaxies can be obtained by appropriately ad- 
justing s°. 

We would like to stress at this point that the 
two panels in Fig. [1] are only shown to illustrate 
that the s° is indeed constrained and to provide 
an idea how the fit quality changes when s° de- 
viates from the best fit s°. The ( x 2 ) 1 ^ 2 vs - s ° 
curves are not used to derive the s° of the cluster 
nor the errors. The estimation of s° and errors 
is done using the histograms of the 5°s obtained 
using the Monte-Carlo simulations shown in Fig. 
[2] and explained below. 

We finally use the best fitting s° of each Monte- 
Carlo run to derive s° for A1689 and estimate the 
errors. In Fig. [2] we show the histograms of s° 
values at which each simulated cluster attains its 
minimum ( x 2 ) 1 ^ 2 ; i- e - the best fitting s° for a 
given Monte-Carlo run. 

As is expected the histograms peak nicely at 
the positions where the mean curves in Fig. Q]also 
have their minima. The strong clustering of the 
histograms between s° = 40 kpc and s° = 90 kpc 
with only a few outliers demonstrates that s° is 
well constrained. The flatter ( x 2 ) 1 / 2 vs. s° curves 
at s° >50 kpc for the NSIE models lead to less well 
defined minima and correspondingly wider distri- 
bution of s° at large s° in the histograms. 

The best fit values of s° for the different de- 
scriptions of the smooth DM component of the 
cluster are shown in Table [U The values given 
are the geometric means of the best fit s° of all 
the simulated clusters. We have used the geomet- 
ric mean to estimate the truncation radius since 
this corresponds to fractional change in mass and 
the ( x 2 ) 1 ^ 2 vs - s ° curves in Fig. [T]are relatively 
symmetric in log(s°) (although not exactly as dis- 
cussed earlier). We also give in Table Q] the esti- 
mated 1- and 2-tr errors of s . The errors have 
been derived from the distribution of the best fit 
s° of the simulated clusters shown in Fig. O For 
this the histograms were interpreted as probabil- 
ity distributions of s° and 1- and 2-er confidence 
intervals were estimated by the regions around the 
mean that contain 68.3 per cent and 95.4 per cent 
of the best fit s° values from the simulation for 
the 1- and 2-cr errors respectively. The asymme- 
try of the distribution becomes evident at higher 
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Fig. 1. — The mean ( x 2 } 1 ^ 2 vs - s ° curve for 1000 simulated clusters for the smooth DM described by both 
NSIE (dotted line, squares) and ENFW (dashed line, triangles), and the combined data from the two smooth 
profiles (solid line, circles). The left panels are calculated for the scaling of the truncation radius with s oc a 
and the right panels for s k <t 2 . The points show the 2-sigma clipped mean for each s°, and the error bars 
show the final sigma of the clipped points. Before clipping, the individual curves were normalized to their 
median ( x 2 ) 1 ^ 2 in order to bring all the curves to a similar ( x 2 ) 1 ^ 2 level for comparison. The median curve 
has been brought back to the level of the mean median ( x 2 ) 1 ^ 2 - The combined curve has been calculated 
from the NSIE and ENFW curves and not from the individual curves for the two models. The minimum 
( x 2 ) 1/2 is obtained at -60-70 kpc. 
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Fig. 2. — The histograms of the s° values at which each simulated cluster attains its minimum. The line 
types used are as above. As is expected, the histograms peak nicely at the positions where the mean curves 
on the top row also have their minima. The strong clustering of the histogram between s° = 40 kpc and 
s° = 90 kpc with only a few outliers demonstrates that s° is strongly constrained. The best fit s° and its 
error has been derived from these histograms. 
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confidence limits as can be seen in the 2-er errors. 

For the scaling law expected theoretically for 
galaxies in clusters ( Merrittl ll983). s = s° X (a/ er°), 



we find s = 64_2g kpc, where the errors given 
are 2-cr errors. For the scaling law used for 
galaxies in the field (s = s° x (a/ c ) 2 ), we find 
s° = 66^26 kpc, the errors are again 2-a errors. 



5. Checks on the Robustness of the Re- 
sults 

We have performed the following checks to con- 
firm that the results presented above are reason- 
able and robust. 



5.1. 



and the Total Mass in Cluster 



Galaxies 



The total mass of a galaxy with a BBS pro- 
file can be easily written in terms of its trun- 
cation ra dius s and vel o city d ispersion a as is 
shown in iBrainerd et all (1996). The total mass 
of the cluster galaxies in A1689 with a scaling 
law for the truncation of the halos of the from 
s = s x (er/er°) Q can be written simply as, 



M tat = 7.3xl0 5 ( f- x x 
kpc 



( <Ji 
\km/s 



(1) 

where sq is the normalization of the scaling law, 
<j° a reference velocity dispersion and <Ji is the 
velocity dispersion of galaxy i. 

In our study we have taken <r° = 220 km/s. 
Note that this er° is only a fiducial value and 
is not related to the L* of the galaxies in the 
cluster. With this cr° and our set of galaxies 
in the cluster, the galaxies have the same total 
mass with the two scaling laws (a=l and a=2) if 
s° a=2 — 0.93 x s a a=l . Note that this relation be- 
tween s° a=1 and s° =2 i s the same for all s° =1 and 
„o 

b a=2- 

The fact that the normalizations s° obtained 
for the two scaling laws are very similar (s° =1 = 
64±^| kpc and s° a=w = 66±i| kpc) for cr°=220 km/s 
provides strong support for the results and our 
analysis. 



5.2. Sensitivity of Cluster Lensing to Ex- 
tensions of Galaxy Halos 

To demonstrate that we are indeed able to mea- 
sure the extension of galaxy DM halos with strong 
lensing we have created clusters with an s° in the 
range [20,80] kpc. For each of these clusters we 
have created a mock set of multiple images that 
are exactly reproduced by the cluster. The mock 
multiple image set is based on the observed mul- 
tiple images so that the cluster setup is as close 
to reality as possible. These clusters with known 
galaxy truncation laws are then analyzed in the 
same way as is done for A1689. 

We find that we are able to recover the input s° 
within a few kpc in all cases. Additionally, both 
the change in the fit quality (A( % 2 } 1/2 -0.2") of 
these new simulated clusters and the distribution 
of the best fit s° is similar to what is observed and 
shown in Figures [1] and [H 

5.3. Effect of the Choice of Multiple Image 
Systems 

We have in addition checked the sensitivity of 
the results to our choice of multiple images. This 
was done by running another 100 Monte-Carlo 
runs with a=l and the smooth DM described by 
NSIE profiles. This time for each Monte-Carlo 
run we selected randomly 20 of the 32 image sys- 
tems to use as constraints for the modeling. The 
multiple images with spectroscopic redshifts were 
always included since they are needed to fix the 
overall mass scale of the cluster. With fewer con- 
straints we obtained essentially the same s° with 
larger spread in the distribution of s° from the 
different runs. The best fit s° obtained with 20 
image systems is s° = 59^ 2 g kpc compared to 
s° = 58^11 kpc with all the image systems, the 
errors are 1-a. 

When only 20 multiple image systems were used 
the absolute ( x 2 ) 1 ^ 2 stayed at the same level as 
with all the 32 image systems. This shows that 
the ( x 2 ) 1 / 2 level is not driven by only a few 
image systems but all image systems contribute 
similarly to the ( x 2 ) 1 / 2 level. The change in 
fit quality between best fit s° and extrema at 
s°=20 kpc and s°=200 kpc with fewer image sys- 
tems is A( x 2 ) 1 / 2 —0.1" showing that also the in- 
dividual s° are less well constrained with fewer 
image systems. That no change in s° is obtained 
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Table 1: Derived s° values and 1- and 2-a errors for s° for the different descriptions of the smooth DM 
component of cluster. The truncation radius s of a galaxy depends on its velocity dispersion a and the 
scaling laws adopted are of the form s = s° x (a/a°) a . The s° values given are the geometric means of 
the individual minima of the simulated clusters. The errors are derived from the distribution of the minima. 
We give both 1- and 2-a errors since the asymmetries of the distributions become more apparent at higher 
confidence limits. The histograms of the minima for the different descriptions of smooth DM are shown in 

Fig. m 
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demonstrates that our results for s° are robust. 

6. Comparison with literature 

The extensions of dark matter halos have been 



Nataraian et a 


. (|l998h. Nataraian et al. 


Gavazzi et al. ( 


20041). and Limousin et al. 



2002), 
2006). 

Strong trunc a tion o f galaxies is found i n 



Nataraian et al. <l 1 998h . Nataraian et al 



and Limousin et al 



(12002ft . 

20061 ) when compared to 
galaxies in the field; the truncation radii of an 
L» galaxy spa n the range 17-55 kpc fo r the 6 clus- 
te rs studied in Nataraian et al. ( 20021 ). The halos 
in iLimousin et al.l ( 20061 ) are truncated more with 
a typical truncation radius below 20 kpc. 

An important difference in the analysis of 



Nataraian et al. | (| 1998. 200 2]) to that of lLimousin et al 
( 20061 ) is that Nataraian et al.l also include strong 



lensing features in the central parts of the clus- 
ters to further constrain the mass profile of the 
selected cluster sample and hence also constrain 
the galaxy halo parameters stronger. This helps 
to better define the shear contribution from the 
cluster galaxies and hence the tr uncation radius 



Anoth er major difference is that ILimousin et al 
( 20061 ) work exclusively with ground based data 
where as Hubble Space Teles cope data are used in 

Natarajan et~aH (|l998l f2002h . 

The large errors in the work of Gavazzi et alj 
are caused by the smoothing scale of 6^=220 kpc/hio 
employed in their analysis which restricts the 



achievable resolution. Although they are not able 
to derive strong limits on the sizes of cluster galax- 
ies, they do find that halos on the periphery of the 
cluster MS0302+17 are more strongly truncated 
than the halos on the central regions of the cluster 
providing thus further confirmation for the tidal 
stripping scenario. 



For the range of a°s for the clusters in lNataraian et al 
our s° is in the range [32,66] kpc for a = 1 
and [16,72] kpc for a = 2 (their s° span 17- 
55 kpc). In Fig. [3] we sho w a comparison between 



our results and tho se of Nataraian et al. (|2002l) 



and ILimousin et al.l (J2006) . For our points we 
also show the scaling of s° with a as dotted and 
dashed lines (a = 1 and a — 2 respectively). The 
lines can be used to convert the s° and errors to 
a a different from 220 km/s, making the com- 
parison between other works easier. The solid line 
shows the s — a pairs for a galaxy with a total mass 
of 5xlO 11 M . The scatter of the points is large, 
though mostly consistent within the large error 
bars. There is some indication that the galaxy 
halos in A1689 are more extended than those in 
most of the other cluster studied. 



Natarajan et alJ (|2002f ) compared their results 



for the density of the cluster at the core radius, 
p(r c ), and the truncation radius of galaxies ob- 
tained in their analysis and fo und results in good 
agreement with iMerrittl (|1983T ). 



S 



-0.5 



40 180km/s U-95 x 10 6 M Q kp C - 3 ) kp °' 

(2) 

Using the res ults for A1689 (0^^^ =1450 km/s 
and r c —77 kpc, lHalkola et alj|2006j ) we get an ex- 
pected truncation radius of 54 kpc, a little smaller 
than the ~65 kpc obtained in this work. This 
(small) difference is in fact also expected since in 
our analysis we measure the truncation radii of the 
galaxies along the line of sight. Some of the galax- 
ies will have large cluster-centric distances despite 
their small projected distances from the center. 
This supports the idea that the galaxy clusters are 
mainly truncated by the tidal fie l d of the global 
potential as assumed by Merrittl ( 1983) and also 



shown in numerical simula tions by iMoore et al 
(ll998l) : lGhigna et all (|2000l) . 



When comparing the results from different 
works it should be noted that weak lensing works 
generally include all the galaxies from the cen- 
ter to the periphe ry of the cluster (although 
iGavazzi et al. 2004 do separate the galaxies in 
radius). This means that the results are averaged 
over the cluster galaxy popul ation out to sev- 
eral Mpc ( Limousin et al. 12006). With our strong 
lensing method we include galaxies only out to a 
projected cluster-centric radius of r ~ 300 kpc. 
The clusters also vary in their central densities 
complicating direct compariso n betw een clusters. 
According to iLimousin et al. ( 20061 ) their clus- 
ter sample (Abell clusters A1763, A1835, A2218, 
A383 and A2390) form a homogeneous set of clus- 
ters and hence the results for these clusters should 
be comparable. 

Comparison to field galaxies is s hown in Fig. |H 



In the f igure we show point s from Braincrd et a 



(199 61), iFischer et alj dloocj). ISmith et al l (l200ll). 



Hoekstral (|2003l ) and lHoekstra et al.l (120041). Adop t- 



ing a° 36 =136 km/s used by Hoekstra et al. ( 20041 ) 



'136 



= 39^17 kpc for a—1 and s 



136 



we obtain 

25^^q kpc for a— 2. Similarly to previous studies 
of cluster galaxies we report a strong truncation of 
galaxy halos in dense cluster environments com- 
pared to galaxy halos in the field. 



± This work 
■ Natarajan 02 
• Limousin 06 




40 60 
s° [kpc] 

Fig. 3. — A comparison of three studies of 
galaxy truncation in the dense cluster environ- 
ments. The red triangles a r e from this work, 
squares fr om [ Nataraj an et alJ (120021) and the cir- 
cles from Limousin et al. (2006). The error bars 



shown are all 1-er. The dotted lines show the scal- 
ing of the best fit values and errors for a = 1 , the 
dashed lines for a = 2. With the solid line we 
show the s — a pairs for a galaxy with a total mass 
of 5xlO n M . 



7. Summary and Conclusions 

In this paper we report the determination of the 
sizes of galaxy dark matter halos in galaxy cluster 
A1689. The strong lensing models for the clus- 
ter are constrained by 107 multiple images and an 
arc in 32 image systems. The strong constraints 
from these images enable us to study not only the 
global mass profile of the clusters but also the ones 
of the cluster galaxies. Assuming well motivated 
scaling laws between the truncation radius of a 
galaxy halo and its central velocity dispersion (as 
obtained with the fundamental plane and Faber- 
Jackson relations) we can study the combined ef- 
fect of the cluster galaxies on the multiple images 
and the ensemble properties of the galaxies. This 
is the first time the sizes of galaxy halos have been 
measured using strong lensing only. 

For a scaling law of the form s ga i = s°x (a ga i /c° 



we find s° 



64 



-67 

-28 



kpc for a = 1 and s° 
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Fig. 4. — A comparison of studies of galaxy trun- 
cation in the field. The red triangles are from this 
work. The dotted lines show the scaling of the best 
fit values and errors for a = 1 , the dashed lines for 
a = 2. The discrepancy is clear between the re- 
sults obtained for galaxies in clusters and those in 
the field. Notice the difference in the horizontal 
scale between this figure and Fig. [31 The solid 
line shows the s — a pairs for a galaxy with a total 
mass of 5xlO n M . 

66^26 kpc for a = 2. Both values are given for a 
fiducial galaxy velocity dispersion of <7°=220 km/s. 
The errors are 2-<r errors to show the clear asym- 
metry of the errors. The s°s are in good agreement 
with previously determined values in several other 
clusters using weak lensing ( Nataraian et al.f l998. 
20021 iLimousin et al.|[2006h . 

Galaxy halos in a cluster can be truncated ei- 
ther by the tida l field of the global c luster potential 
or harassment ( Moore et al. 19961 1998) by other 
cluster galaxies that strip the halos of galaxies in 
the central regions of cluster. Mergers of clus- 
ter galaxies on the o ther hand are extremely rare 



(jGhigna et alj 119981 ). Once the cluster has been 



formed the principal mechanism for truncation is 
the tidal strippi ng of galaxy halos by the global 
cluster potential ( Ghigna et alj200Cl) . This is sup- 
ported by the correlation between the density of 
the cluster at the core r adius and the truncatio n 
radii of galaxies shown in lNatarajan et al ] (|2002h . 



stripping scenario. 
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